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摘 要 : 利用 较 高 空间 分 辨 率 的 Landsat 卫星 数据 ,以 覃 志 豪 单 窗 算法 反 演 黄河 源 区 1990 一 2018 年 
(间隔 3~Sa) 的 夏季 (7 月 或 8 月 ) 地 表 温 度 ,并 对 反 演 结果 进行 相关 性 分 析 。 结 果 表 明 :(1) 黄河 源 
区 夏季 地 表 温 度 平均 值 与 MODIS 地表 温度 产品 的 平均 值 相 关系 数 不 高 , 约 为 0.5, 主 要 原因 是 
Landsat 数据 与 MODIS 产 品 的 空间 分 辩 率 不 一 致 .(2) 水 体 的 地 表 温 度 基 本 保持 不 变 ,其 他 地 物 的 
地 表 温 度 均 呈现 上 升 趋势 ,冰川 的 地 表 温 度 升 高 最 快 ,(3) 黄河 源 区 各 年 份 地 表 温 度 与 高 程 存在 明 
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于 -0.5 的 区 域 约 占 整个 源 区 面积 的 70%。 


显 的 负 相 关 关 系 ,相关 系数 平均 值 为 -0.65。.(4) 空间 上 ,黄河 源 区 地 表 温 度 值 与 土壤 湿度 值 的 空间 
分 布 存 在 负 相 关 关 系 ,(5) 在 气候 因素 中 ,降水 与 地 表 温 度 存 在 明显 的 负 相 关 关系 ,其 相关 系数 小 
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黄河 源 区 作为 黄河 的 发 源 地 ,是 中 华 水 塔 的 重 
要 组 成 部 分 。 源 区 生态 环境 敏感 .脆弱 , 近 几 十 年 
来 , 源 区 多 年 冻 土 已 发 生 大 面积 退化 并 日 趋 严重 ， 
植被 范围 逐渐 减少 ,草场 退化 水 土 流失 严重 等 生 
态 环境 问题 日 趋 严 重 。 地 表 温 度 是 地 表 大气 及 太 
阳 之 间 的 热量 交换 ,不 少 学 者 将 其 视 为 监测 地 球 环 
境 动 态 变 化 的 重要 指标 之 一 。 有 学 者 将 地 表 温 度 
应 用 于 灾害 监测 ,如 利用 热 红 外 异常 来 研究 地 表 温 
度 与 地 震 发 生 的 关系 "。 研 究 黄 河源 区 地 表 温 度 是 
深入 研究 三 江 源 地 区 能 量 循环 .环境 保护 的 基础 ， 
可 为 黄河 源 区 的 生态 保护 及 三 江 源 气候 变化 与 生 
态 响 应 提供 科学 依据 。 同 时 ,精确 反 演 地 表 温度 对 
于 全 球 气候 变化 的 研究 地球 环境 灾害 监测 及 自然 
资源 管理 有 着 极其 重要 的 意义 ”。 

单 通道 算法 是 利用 卫星 传 感 顺 上 的 一 个 热 红 
外 通道 来 获得 地 表 的 辐射 能 量 ,并 利用 大 气 的 温度 
ERR .湿度 廓 线 和 压力 廓 线 ,建立 大 气 辐射 传输 方 
程 ,以 此 为 基础 来 消除 大 气 以 及 地 表 比 辐射 率 对 地 
表 温 度 的 影响 ,进一步 通过 普 朗 克 函 数 反 演 得 到 地 
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表 温 度 ”。 覃 志 豪 等 "利用 Landsat TM 数据 的 第 6 
波段 ,推出 了 只 需 3 个 参数 反 演 地 表 温 度 的 算法 ,其 
具有 和 较 高 反 演 精度 。 毛 克 彪 等 5 首次 利用 重 志 察 提 
出 的 基于 Landsat TM 数据 的 单 窗 反 演 算法 ,对 MO- 
DIS 卫 星 的 第 31、32 波 段 建 立 了 辐射 率 与 温度 的 线 
性 经 验 关 系 式 。 毛 克 彪 等 “首次 推导 出 了 适用 于 
ASTER 传感器 的 单 窗 算 法 ,该 算法 更 适合 于 小 区 域 
的 温度 空间 差异 分 析 。Landsat 卫星 具有 空间 分 辩 
率 高 .观测 时 间 序 列 长 的 优势 ,已 经 有 接近 40 a 的 观 
测 数据 。 其 次 ,Landsat 是 一 种 多 光谱 数据 ,光谱 信 
息 丰 富 且 存储 数据 多 ,其 热 红外 波段 对 地 物 的 热 信 
息 比 较 敏 感 。 因 此 ,本 文选 用 Landsat 数 据 反 演 黄 河 
源 区 夏季 的 地 表 温 度 。 


1 数据 与 方法 


1.1 研究 区 概况 
黄河 源 区 (33°~35°N,95°~102°E) 位 于 青藏 高 
原 东 部 .青海 省 南部 (图 1) ,南北 长 约 171 km, 东 西 
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图 例 
县 城 驻 地 0 100 km 


图 1 研究 区 示意 图 
Fig. 1 Schematic diagram of the study area 


长 约 232 km ,面积 约 14.75 km ,平均 海拔 在 4000 m 
左右 ,最 高 约 6200 m。 研 究 区 属 高 原 大 陆 性 气候 ， 
气温 相对 较 低 ,早晚 温差 大 ,夏季 多 年 平均 气温 为 
3.67 CC。 源 区 蒸 散 发 相对 强烈 ,年 均 莱 发 量 达 到 
1200~2000 mm,2001 一 2009 年 黄河 源 年 蒸 散 量 增加 
明显 。 
1.2 研究 数据 
本 研究 选取 1990 1995 , 2000, 2005 , 2010 年 7 
月 或 8 月 Landsat5 的 TM 数据 ,选取 2015、2018 年 7 
月 或 8 月 Landsat8 的 OLITIRS 数 据 ,由 于 Landsat 扫 
描 面积 相对 较 小 ,同一 月 份 很 难 拼 出 完整 覆盖 研究 
区 的 有 影像 ,需要 多 幅 影 像 馈 般 而 成 ,选取 时 云 量 低 
F 10% , 则 向 相 邻 (前 后 ) 年 度 取 缺 少 部 分 的 影像 。 
本 文采 用 空间 分 辩 率 为 30 m 的 SRTM-DEM ATEX 
据 ,来 自 中 科 院 地 理 空间 数据 云 网 站 (http://www.gs- 
cloud.cn) ;气象 数据 来 自 中 国 气 象 数 据 网 (http://da- 
ta.cma.cn) ;土壤 湿度 数据 来 自 GLDAS 全 球 陆 面 数 
据 同化 系统 (https://disc.gsfc.nasa.gov/datasets ) o 
1.3 研究 方法 
13.1 反 演 方法 ， 覃 志 豪 针对 Landsat 5 数据 的 热 红 
外 波段 特点 ,根据 地 表 热 辐射 传导 方程 推导 出 了 精 
度 较 高 的 地 表 温 度 反 演算 法 ,主要 是 将 大 气 和 地 
表 的 作用 直接 反映 在 公式 中 ,公式 如 下 4 
T={al -C-D)+[b(1-C-D)+C+DIf,,,, -DT JC (1) 
C=te (2) 
D=(1 -D1 +(1 -e)r] (3) 
式 中 ;7 为 实际 地 表 温 度 (K);a、b 为 回归 系数 (a= 
-37.3556,b=0.4586) ; T na 为 传感器 亮度 温度 (K); 
7 为 大 气 平均 作用 温度 (K); t 为 大 气 透 过 率 ; e 为 
地 表 比 辐射 率 ;C.D 为 关于 7 与 的 孔 数 。 该 方法 
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适用 于 只 有 一 个 热 红外 波段 的 遥感 数据 , 且 仅 需 
地 表 比 辐射 率 、 大 气 透 过 率 、 大 气 平均 作用 温度 3 个 
参数 即 可 实现 地 表 温 度 的 反 演 。 
国内 外 学 者 针对 Landsat8 TIRS 的 第 10 波 段 , 对 
单 窗 算 法 公式 中 的 ab 2 个 参数 进行 拟 合 ,得 到 
Landsat8 TIRS 波 段 10 的 地 表 温 度 反 演算 法 公式 中 : 
T,= {a(l - Cio = Dio) +[bio(l -= Cio = Da) + 
Cio t Diol sso = Dial Cio (4) 
式 中 : T, 为 Landsat8 TIRS 波 段 10 反 演 得 到 的 地 表 
温度 (K); TETA 为 Landsat8 TIRS Y Be 10 WERA 
亮度 温度 (K); ao 、bw 为 回归 系数 ( 表 1);Co、Dn 
为 关于 中 间 变 量 to 与 21 的 函数 ,根据 以 下 公式 得 
到 : 
CoE Tow Di=(L-znol+roL-eo ©) 


#1 Landsat8 TIRS 波 段 10 的 回归 系数 取 值 ”” 
Tab.1 Regression coefficient value of Landsat8 
TIRS band 10" 


Yak BE/SC aio bio R 
20~70 -70.1775 0.4581 0.9997 
0~50 62.7182 0.4339 0.9996 

-20~-30 -55.4276 0.4086 0.9996 


注 :aw .bw 为 Landsat8 TIRS 波 段 10 的 回归 系数 ;RR 为 决定 系数 。 


1.3.2 分 析 方 法 ”本文 分 析 利 用 一 元 线性 回归 趋势 
分 析 方 法 来 计算 每 个 像 元 在 1990 一 2018 年 (间隔 3~ 
5 a) 的 变化 趋势 ,分 析 黄 河源 区 夏季 地 表 温 度 的 变 
化 情况 ,计算 公式 为 


nx Dix QinLST,- Hid ix QinLST, 
=A i=l i=l 


QinLST ppe = (6) 


式 中 : QinLST,,, 为 每 个 像 元 在 这 7 a 变化 趋势 的 斜 
率 ;i 为 年 序号 (j=1, 2, 3,…,n); QinLST, 为 第 i 年 夏 
FWE as Se IK Hl Ze lid EEL. QinLST,,.>0 说 明 
地 表 温 度 在 上 升 , QinLST,, .<0 ,说 明 地 表 温 度 值 
在 下 降 。 

数学 家 皮尔 逊 基 于 极 大 似 然 法 提出 变量 之 间 
的 相关 系数 (7) ,被 用 来 衡量 2 个 变量 之 间 线 性 关系 
的 强 弱 。r 的 正 负 代表 着 变量 间 的 相关 性 质 , 若 其 " 
越 接 近 于 1, 则 说 明 两 组 数据 之 间 的 正 相 关 性 越 
强 。 本 文通 过 计算 皮尔 逊 相关 系数 来 反映 地 表 温 
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度 与 其 影响 因素 之 间 的 关系 ,计算 公式 如 下 
Di - 9) x0,-7) 
= = - (7) 
Se- So 

式 中 :i 为 第 几 年 ;n 为 时 间 序 列 长 度 ;xi yi 为 第 i 年 的 
样本 值 ; 去 .7 了 分别 为 样本 值 x、y 的 平均 值 。 
1.3.3 反 演 参数 的 确定 蕴 志 坚 地 表 温度 单 窗 反 演 
算法 需要 确定 3 个 物理 参数 :(1) 地 表 比 辐射 率 , 指 
物体 与 黑体 在 同 温度 、 同 波长 下 辐射 出 射 度 的 比 
值 * 中 。 在 确定 传 感 融 测 中 范围 和 像素 尺寸 时 ,地 
表 比 发 射 率 的 大 小 主要 取决 于 地 表 材 料 的 组 成 . 形 
状 和 结构 "中 。 研 究 区 的 地 表 比 辐射 率 根据 其 像 元 
中 的 土地 利用 类 型 及 植被 指数 计算 "“"”, 本 研究 区 
将 地 表 划 分 为 水 体 、 植 被 .自然 地 表 ( 裸 地 + 建筑 用 
地 + 冰川 ).(2) 大 气 透 过 率 , 指 电磁 波 在 大 气 中 传播 
时 ,经 大 气 衰减 后 的 电磁 辐射 通 量 与 人 射 时 电磁 辐 
射 通 量 的 比值 ,其 大 小 主要 由 大 气 状况 入 射 点 高 
度 等 因素 决定 。 本 文 利用 草 志 肾 等 提出 结合 水 汽 含 
量 动态 变化 的 大 气 透 过 率 方 程 进行 估算 (3) 大 
气 平 均 作 用 温度 ,主要 受 大 气 分 子 的 活跃 状态 .大气 
剖面 温度 及 大 气 成 分 的 运动 特征 等 影响 所 。 本 文 将 
根据 大 气 平均 气温 与 地 表 高 程 之 间 的 关系 ,通过 大 
气 模拟 ,得 出 在 不 同 大 气 模式 及 不 同 大 气 剖 面 的 背 
景 下 ,估算 大 气 平均 作用 温度 ”站 。 其 中 地 表 比 辆 
射 率 的 计算 需要 用 Landsat 多 光谱 数据 求 得 植被 指 
AX NDVI, 再 利用 像 元 二 分 模型 计算 得 到 植被 覆盖 
率 , 结 合 不 同 地 物 类 型 的 地 表 求 得 整个 源 区 地 表 比 
辐射 率 。 大 气 透 过 率 主要 受 大 气 水 汽 含量 的 影响 ， 
借助 经 验 公 式 求 得 大 气 透 过 率 的 值 。 根 据 不 同 大 
气 模式 的 经 验 公 式 估算 得 到 大 气 平均 作用 温度 。 
1.3.4 地 表 温 度 反 演 ”原始 数据 经 过 一 系列 预 处 理 
后 ,需要 进行 辐射 亮 温 的 计算 ,利用 Planck 公式 进 
行 亮 度 温 度 值 的 转化 , 即 可 获得 地 表 温 度 ,公式 如 
T; 


K 


了 = 一 一 一 一 一 (8) 
sensor l 1 $ K, 
a 了 
式 中 :Tow 为 传 感 带 亮 度 温度 (KK) ;Ze 为 经 过 预 处 


理 之 后 的 辐射 亮度 ( wm2sr bm ); K,、K, 为 常 
量 ,可 从 影像 数据 文件 中 获得 ( 表 2)。 
根据 公式 (1) 、(4) (6) 即 可 反 演 得 到 地 表 


表 2 辐射 常数 K 


Tab.2 Radiation constant K 


卫星 及 波段 K, K: 
Landsat5 波段 6 607.76 1260.56 
Landsat8 波段 10 774.89 1321.08 


2 结果 与 分 析 


2.1 反 演 结果 的 验证 

将 硬 志 豪 反 演 地 表 温 度 (简称 反 演 地 表 温 度 ) 
与 MODIS 地 表 温 度 进 行 对 比 验 证 。 由 于 MODIS 数 
据 仅 在 2000 年 以 后 才 有 数据 产品 ,因此 选取 2000、 
2005.2010.2015 .2018 R Æ MODIS 地表 温度 数据 
进行 验证 (图 2)。 

总 体 上 , 反 演 地 表 温 度 的 平均 值 高 于 MODIS 地 
表 温 度 的 平均 值 , 这 是 由 于 相对 于 其 他 热 红 外 传 感 
an , Landsat 数据 具有 非常 高 的 空间 分 辨 率 , 而 MO- 
DIS 数据 分 辩 率 较 低 ,其 混合 像 元 包含 的 地 物 信 息 
变化 相对 更 大 ,会 拉 低 整 体 的 地 表 温 度 值 ,另外 黄 
河源 区 扎 陵 湖 及 哪 陵 湖 的 面积 较 大 ,进一步 影响 了 
黄河 源 区 夏季 地 表 温 度 的 平均 值 。5 个 年 份 的 地 表 
温度 平均 值 在 2015 年 发 生 突变 ,MODIS 地 表 温 度 
平均 值 在 2015 年 之 前 是 先 下 降 而 后 再 上 升 ,而 反 
演 地 表 温度 在 2015 年 之 前 呈现 先 上 升 再 下 降 而 后 
再 上 升 的 趋势 ,在 2015 年 之 后 ,两 种 数据 地 表 温 度 
的 平均 值 呈 下 降 趋 势 。 在 年 际 尺度 上 ,两 种 地 表 温 
度 平 均值 有 所 差异 ,因此 有 必要 对 两 种 地 表 温 度数 
据 进 行 相关 性 分 析 。 


D 反 演 地 表 温度 ” 国 MODIS 地 表 温 度 


地 表 温度 /°C 


2000 2005 2010 2015 2018 
年 份 


图 2 反 演 地 表 温 度 与 MODIS 地 表 温度 对 比 
Fig.2 Comparison between retrieved surface temperature 


and MODIS surface temperature 
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2000—2018 年 反 演 地 表 温 度 与 MODIS 地 表 温 
度 呈 正 相 关 , 其 中 2015 年 两 种 地 表 温度 数据 的 相关 
性 最 高 ,r 为 0.59,2018 年 两 种 地 表 温 度 的 相关 性 最 
低 , 其 余 各 年 份 - 值 不 是 很 高 ,在 0.43~0.53 之 间 ( 图 
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3) ,主要 原因 为 两 种 地 表 温 度数 据 的 空间 分 辨 率 不 
一 致 ; 反 演 的 Landsat 数 据 是 由 不 同时 间 的 影响 锐 散 
而 成 ,存在 一 定 误 差 ;少量 云 的 存在 会 影响 反 演 的 


60 (a) 2000 年 60 r (b) 2005 年 60 F (c) 20104 
50 at Be 50 50+ 3 
o 40 o 40 2 40r 
网 30 = 30 = 30+ 
aK 
p 20 a 20 : 0.51 8 20+ 
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图 3 反 演 地 表 温度 与 MODIS 地 表 温 度 的 相关 性 


Fig. 3 Correlation between retrieved surface temperature and MODIS surface temperature 


2.2 地 表 温 度 的 时 空 分 布 特征 

利用 覃 志 豪 单 窗 反 演 算法 反 演 得 到 黄河 源 区 
夏季 地 表 温 度 空间 分 布 ( 图 4a~g)。 整 体 上 ,黄河 源 
区 北部 地 表 温 度 最 高 ,可 达 50 % ,这 是 因为 北部 地 
区 海拔 较 低 , 且 耕 地 及 建筑 用 地 较 多 ; 黄河 源 中 部 
地 区 地 表 温 度 呈 现 蓝 色 的 区 域 为 冰川 ,其 温度 最 
低 , 约 为 -24 % ,其 他 蓝 色 的 区 域 有 部 分 是 云层 覆 
盖 , 导 致 地 表 温 度 较 低 呈 现 蓝 色 。 

局 部 来 看 ,黄河 源 区 西部 的 两 大 湖泊 , 扎 陵 湖 
及 哪 陵 湖 的 表面 温度 相对 较 低 ;1995 年 黄河 源 区 南 
部 .2005 年 及 2010 年 黄河 源 区 西南 部 地 表 温 度 有 负 
值 出 现 (图 4b .4~e) ,由 于 零散 的 云 覆 盖 所 导致 。 

根据 线性 回归 公式 (6) 计 算出 7 个 年 份 地 表 温 
度 的 变化 斜率 值 ( 图 4h) ,对 地表 温度 年 平均 值 序列 


变化 进行 趋 热 分析, 从 图 中 可 以 看 出 ,1990 一 
2018 年 ,黄河 源 区 北部 及 中 部 边缘 山区 地 表 温 度 呈 
现 上 升 趋势 ,黄河 源 西 部 及 东部 地 区 地 表 温 度 呈 现 
下 降 趋 势 。 
2.3 地 表 温度 影响 因素 分 析 

地 表 温 度 主 要 受 高 程 .土壤 湿度 .气候 、 植 被 状 
况 、 用 地 类 型 以 及 人 类 活动 等 因素 的 共同 影响 。 黄 
河源 区 建筑 用 地 少 ,人 类 活动 范围 较 小 ,因此 黄河 
源 区 地 表 温 度 主要 受 自 然 因 素 影 响 。 降 水 和 气温 
通过 水 热 交换 影响 了 每 个 像 元 地 表 温 度 的 大 小 ;在 
众多 影响 因素 中 ,地表 高 程 与 地 表 温 度 的 相关 性 最 
高 ;土壤 湿度 也 影响 了 地 表 温 度 的 空间 分 布 ,植被 
及 土地 利用 类 型 对 地 表 温 度 的 影响 主要 是 通过 地 
表 比 辐射 率 体现 的 所 。 因 此 ,探讨 高 程 .土壤 湿度 、 
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图 4 反 演 地 表 温 度 的 空间 分 布 及 年 际 变化 趋势 


Fig.4 Spatial distribution and interannual variation trend of retrieved surface temperature 


气候 ,植被 状况 和 土地 类 型 等 因素 对 黄河 源 区 夏季 
地 表 温 度 时 空 分 布 的 影响 十 分 必要 。 

2.3.1 高 程 根据 研究 区 高 程 数 据 ,选取 高 程 起 伏 
较 大 ,纬度 变化 小 .地 物 类 型 变化 不 大 的 区 域 ,以 保 
证 在 高 程 与 地 表 温 度 的 相关 性 分 析 中 ,高 程 是 最 主 
要 的 影响 因素 。 提 取 其 7 个 年 份 的 高 程 以 及 地 表 温 


土壤 湿度 的 相关 关系 。 但 是 黄河 源 区 面积 大 ,可 以 
在 空间 上 进行 对 比 ,可 在 一 定 程度 上 反映 土壤 湿度 
与 地 表 温 度 的 变化 情况 。 黄 河源 区 土壤 湿度 的 高 
值 区 分 布 在 研究 区 的 南部 及 东南 部 ,而 低 值 区 域 分 
布 在 东北 部 及 西部 (图 6a)。 空 间 上 ,黄河 源 区 东部 
及 西部 地 表 温 度 较 高 ,其 土壤 湿度 低 ;黄河 源 区 东 


度 并 计算 其 > 值 ( 图 5) ,各 年 份 地 表 温 度 与 高 程 存在 
明显 的 负 相 关 关 系 ,r 值 平均 为 -0.65,2000、2015、 
2018 年 均 达 0.70 以 上 (图 5c、f~g)。 高 程 平均 每 增 
高 100 m, 地 表 温 度 下 降 约 0.925 Co 

2.3.2 土壤 湿度 ”选用 GLDAS 土壤 湿度 数据 产品 
(https://disc.gsfc.nasa.gov ) ,该 数据 的 空间 分 辨 率 较 
为 粗糙 , 仅 为 0.25" ,无 法 在 空间 上 反映 地 表 温 度 与 


南部 地 表 温 度 相 对 较 低 ,其 土壤 湿度 高 , 某 种 程度 
上 可 以 说 明 地 表 温 度 与 土壤 湿度 在 空间 分 布 上 存 
在 负 相 关 关 系 (图 6)。 

2.3.3 气温 将 气温 栅 格 数据 重新 降 尺 度 至 30 m 
空间 分 辨 率 , 并 投影 至 与 地 表 温 度 空间 分 布 图 一 至 
的 坐标 系 , 计 算出 黄河 源 区 夏季 地 表 温 度 与 夏季 气 
温 的 相关 性 空间 分 布 (图 7) ,整个 黄河 源 区 地 表 温 
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Fig. 5 Correlation between surface temperature and elevation 
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图 6 2015 年 土壤 湿度 与 地 表 温 度 的 空间 分 布 


Fig.6 Spatial distributions of soil moisture and surface temperature in 2015 


度 与 气温 呈正 相关 的 部 分 约 占 55%,r 平 均值 约 为 
0.38, 主要 分 布 在 黄河 源 区 的 西部 .北部 和 东南 部 山 
区 。 源 区 西部 及 东部 存在 明显 的 分 带 现 象 , 是 由 于 
不 同时 相 的 影像 镶 谍 而 成 。 

2.3.4 降水 量 利用 ENVI 软 件 计 算得 出 黄河 源 区 夏 
季 地 表 温 度 与 夏季 降水 量 的 相关 性 空间 分 布 (图 8)， 
进一步 统计 顶 格 ,其 中 约 70% 的 像 元 其 > 为 负 值 , 其 
中 r+ 小 于 -0.4 的 像 元 约 占 20% ,说 明 随 着 该 地 区 降 
水 量 的 增加 ,相应 的 黄河 源 区 地 表 温 度 降 低 , 呈 负 
相关 关系 的 像 元 主要 分 布 在 黄河 源 区 西部 两 大 湖泊 
附近 ,及 黄河 源 中 部 山区 、 南 部 东部 边缘 地 区 。 
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图 7 地 表 温 度 与 气温 的 相关 性 
Fig.7 Correlation between surface temperature and 


air temperature 
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图 8 地 表 温度 与 降水 量 的 相关 性 
Fig. 8 Correlation between surface temperature and 
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2.3.5 用 地 类 型 ”提取 不 同 地 物 类 型 的 地 表 温 度 
(图 9) ,不 难看 出 ,1990 一 2018 年 水 体 的 地 表 温 度 基 
本 不 变 ,其 他 地 物 的 地 表 温 度 均 呈现 上 升 的 趋势 ; 
冰川 与 襟 地 的 地 表 温 度 升 高 最 快 ,冰川 从 1990 年 平 
均 地 表 温 度 -7.0 CE Ft Bl) 2018 年 平均 地 表 温 度 
3.0 % ,增幅 达 10.0 °C ; 裸 地 由 1990 年 平均 地 表 温 度 
23.7% 上升 到 2018 年 的 34.5 C ,增幅 达 10.8 °C; #F 
地 植被 的 地 表 温 度 呈 现 上 升 趋势 ,但 增 速 较 缓 ,分 
别 由 1990 年 平均 地 表 温 度 27.2 C 21.2 % 上升 到 
2018 4E AY) 34.1 C .27.3 % ,增幅 分 别 达 6.9 % .6.1 Co 


国 植被 E 水 体 冰川 
2005 2010 2015 2018 
年 份 


图 9 用 地 类 型 与 地 表 温度 年 际 变化 


Fig.9 Interannual variation of land use type and surface temperature 


3 结论 


i 


本 文选 用 7a 的 Landsat TÆR He , Fl) HY aR BE 
单 窗 算法 反 演 黄河 源 区 1990 一 2018 年 (间隔 3~5 a) 
夏季 的 地 表 温 度 ,通过 一 元 线性 回归 及 相关 性 分 析 
方法 ,分 析 得 到 黄河 源 区 地 表 温 度 的 时 空 分 布 特征 
和 变化 趋势 及 影响 因素 ,结论 如 下 : 

(1) 1990 一 2018 年 研究 区 北部 及 中 部 边缘 山区 
地 表 温 度 呈 现 上 升 趋势 ,研究 区 西部 及 东部 地 区 地 
表 温 度 呈 现下 降 趋 势 。 

(2) 1990—2018 年 水 体 的 地 表 温 度 基 本 不 变 ， 
其 他 地 物 的 地 表 温 度 均 呈现 上 升 的 趋势 ;冰川 与 裸 
土 的 地 表 温 度 升 高 最 快 ,增幅 分 别 达 10.0 %C、 
10.8 ;耕地 及 植被 的 地 表 温 度 值 增长 速度 较为 组 
慢 , 由 于 黄河 源 区 植被 以 草地 为 主 ,地 表 温 度 对 植 
被 覆盖 变化 的 响应 存在 一 定 的 正 相 关 关 系 。 

(3) 研究 区 高 程 与 地 表 温 度 + 为 -0.65, 即 高 程 


平均 每 升 高 100 m, 地 表 温 度 约 下 降 0.925 °C ; 黄河 
源 区 东部 及 西部 地 表 温 度 较 高 ,相应 的 其 土壤 湿度 
较 低 ;黄河 源 区 东南 部 地 表 温 度 相 对 较 低 ,其 土壤 
湿度 较 高 。 

(4) 气温 的 分 布 对 地 表 温 度 的 时 空 分 布 也 有 
一 定 影响 ;降水 量 与 地 表 温 度 呈 负 相 关 关 系 ,黄河 
源 区 约 70% 的 像 元 其 > 为 负 值 ,说 明 随 着 该 地 区 降 
水 量 的 增加 ,相应 区 域 的 地 表 温 度 降低 , 呈 负 相关 
关系 的 像 元 主要 分 布 在 黄河 源 区 西部 两 大 湖泊 附 
近 , 及 黄河 源 中 部 山区 .南部 东部 边缘 地 区 。 
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Change of surface temperature in the source area of the Yellow River in summer 
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Beijing, School of Water Resources and Environment, Beijing 100083, China) 


Abstract: Using Landsat satellite data with high spatial resolution, the Qin Zhihao mono window algorithm was 
used to retrieve the summer (July or August) surface temperature in the source area of the Yellow River from 
1990 to 2018 (with an interval of 3-5 years), and correlation analysis was conducted on the retrieval results. The 
results show that: (1) The correlation coefficient between the average retrieved surface temperature and the aver- 
age MODIS surface temperature product is not high, about 0.5, mainly due to the inconsistent spatial resolution 
between Landsat data and MODIS products. (2) The surface temperature of water remains basically unchanged, 
while the surface temperature of other features shows an upward trend, with glaciers experiencing the fastest in- 
crease in surface temperature. (3) There is a significant negative correlation between surface temperature and ele- 
vation in the source area of the Yellow River in various years, with an average correlation coefficient of - 0.65. 
(4) In terms of space, there is a negative correlation between the spatial distribution of surface temperature values 
and soil moisture values in the source area of the Yellow River. (5) Among climate factors, there is a significant 
negative correlation between precipitation and surface temperature, with regions with a correlation coefficient 
less than -0.3 accounting for 70% of the entire source area. 

Key words: the source area of the Yellow River; surface temperature; remote sensing retrieval; changing char- 


acteristics; controlling factors 


